EXPONENTIALLY LONG TIME STABILITY NEAR AN 
EQUILIBRIUM POINT FOR NON LINEARIZABLE ANALYTIC 

VECTOR FIELDS. 



TIMOTEO CARLETTI 



Abstract. We study the orbit behavior of a germ of an analytic vector field 
of (C",0), n > 2. We prove that if its linear part is semisimple, non-resonant 
and verifies a Bruno-like condition, then the origin is effectively stable: stable 
for finite but exponentially long times. 



1. Introduction 

Let us consider the germ of analytic vector field, Xp — J2i<j< n ( z )~{Tt~> °^ 
(C n jO) n > 2, whose components (Fj)i<j< n are analytic functions vanishing at 

o e c™. 

Let us consider the associated Ordinary Differential Equation: 

dz , , 

(1-1) Tt =F(z); 

under the above assumptions z(t;0) = for all t is an equilibrium solution . We 
are interested in studying the stability of orbits of Xp in a neighborhood of this 
equilibrium point. 

We use the standard definition of stability (see |Mop for an equilibrium solution: 
z = it is stable is the past and in the future if for any neighborhood U of there 
exists a neighborhood V, containing the origin, s.t. z(0; zq) G V implies z(t; zq) G U 
for all t ER. 

In a coordinates system centered at the equilibrium point the j'-th component of 
the vector field will take the form: Fj(z) = (Az)j + fj{z), with A & n x n complex 
matrix and fj analytic function such that fj(0) = Dfj(Q) = 0, for all 1 < j < n. 
Following the idea of Poincare to study the orbit of in a neighborhood of the 
origin, we will try to find an analytic change of coordinates, through an analytic 
diffcomorphisms z H{z) = w the linearization, s.t. in the new coordinates the 
vector field Xp is conjugate to its linear part, Xa — ~^2,{Az)j-^-: H* XpH~ x = Xa- 
Hence equation 1)1. 1|) rewrites: 

t * dw 

(1.2) -=A W . 
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^Here and throughout the paper by z(t; zq) we mean the solution at time t of 11.11 s.t. 2(0; zq) = 
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This change of coordinates must solve: 

(1.3) AH(z) = DH{z) ■ (Az + f{z)) , 

and it is unique by assuming DH(0) = I. 

Clearly if the linear system 11.2fl is stable and (|l.lfl is analytically linearizable, 
then also the latter is stable. It is a remarkable result that this condition is also 
necessary, as the following Theorem states: 

Theorem 1.1 (Caratheodory-Cartan 1932). Necessary and sufficient condition for 
the stability of the solution z = of ()1.1J) for all real t is that: 

(1) A is diagonalizable with purely imaginary eigenvalues; 

(2) there exists an holomorphic function z — K(w) = w + 0(\w\ 2 ), w G C", 
which brings (11.11) into the linear system: 

dw 

— = Aw . 

dt 

So let us assume A to verify hypothesis of Thcorcm ll.il let (u>j)x<j< n C M and 
A = diag(iuji , . . . , iuj n )- 

Then A belongs to the Siegel domain 2 : the origin is contained in the convex 
hull of the set of eigenvalues plotted as points in the complex plane (e segment in 
this case). This is the harder situation w.r.t. to the complementary case, Poincare 
domain, because small divisors are involved: the existence of an analytic lineariza- 
tion is strictly related to the arithmetic property of approximation of the vector 
u> = (u>i, . . . , uj n ), with vectors with integer entries. 

The first step is to assume A to be non-resonant: a ■ uj ^ uij, for all a G N™ s.t. 
|a| = ai + • • • + a n > 2 and for all j G {1, . . . , n}. This ensures the existence of a 
formal change of variable which linearizes (|l.l|l . 

In |Br| author introduced the, today called, Bruno condition 3 to characterize 
the rate of approximation of vectors of R n with vectors of Z™ : 

nA , ^ logn-^fc + i) ^ 

(1.4) }^ Y k + °°' 

fe>0 

where for all positive integers k, fi(fc) = min{|a • uj — u)j\ : j G {1, . . . , n},a G 
Z™, < \a\ < 2 k }. From the result of Bruno |Brj follows that if A is non-resonant, 
diagonal with purely imaginary eigenvalues and verifies the above condition, then 
there exists an analytic linearization which brings l|l.l|l into (|1.2|) . 

Let us make a step backward and consider the following problem |Ca2| . Let 
Ai C A2 C C" [[zi, . . . , z n }] be two classes of formal power series closed w.r.t. to 
derivation and composition. Let / G A\ s.t. / = J2\ a \>2 fa z °i ^ A G GL(n,C) 



2 According to the classification of |Br| this case is Poincare domain l.d, but we prefer con- 
sider it as a Siegel case because the obstructions to the linearizability are very similar to those 
encountered in the Siegel domain. 

3 The Bruno condition can be rewritten using a general increasing sequence of integer numbers, 
(Pfc)fe- In IBrl pag. 222, author proved that 11.41 is equivalent to: 

y logn-Wx) <+QO; 

where f2(p) = min{|a • w — \ : j £ {1, . . . , n}, a £ Z n , < \a\ < p, }. 
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and consider the following (formal) ODE: 

(1.5) ^ = Az + f(z). 

We say that l|1.5(l is linearizable in A2 if there exists H G A%, normalized with 
H = z + C(|z|l Q l), \a\ > 2, s.t. formally we have: 

w = H(z) and — — = Aw . 

If both A\ and A2 coincide with the ring of formal power series we already 
know that generically the problem has solution if and only if A is non-resonant, 
which will be assumed from now. In the other cases of general algebras of formal 
power series, new arithmetical conditions on A have to be imposed if we are in 
the Siegel domain. This case has been considered in details in |Ca2j section 5, to 
which we refer for all details 4 . There author proved that the Bruno condition is 
still sufficient to linearize whenever A\ = A2, otherwise new Bruno-like conditions 
are introduced, weaker than the original Bruno condition. 

An interesting case is when A\ is the ring of convergent power series in some 
neighborhood of the origin, and A2 is the algebra of Gevrey-s, s > 0, formal power 
series. Namely we are considering the Gevrey linearization of analytic vector fields. 

Let F = J2 f a z a , (/ a )a£N n C C™ be a formal power series, then we say that it 
is Gevrey-s jBal IRaj . s > 0, if there exist two positive constants C\, C2 such that: 



(1.6) |/ | <CxC 2 - sH (|a|O s VaeN". 

We denote the class of formal vector valued power series Gevrey-s by C s . It is 
closed w.r.t. derivation and composition. 

In the Gevrey-s case the arithmetical condition introduced in |Oa2| . called 
Bruno-s condition, s > 0, for short B s , reads: 

(1.7) hmsup ( 2 V lQS "" 1(Pm+l) - slog H ) < +00 , 

|a|-+oo \^ m=0 Pm J 

for some increasing sequence of positive integer (j)k)k and k(o) is defined by p K ( a ) < 

|"| < P«(a)+1- 



Remark 1.2. This definition recall the classical one of Bruno |Br) . where first 
one suppose the existence of a strictly increasing sequence of positive integer such 
that (|1.7I) holds, then one can prove (see |Br| §IV page 222) that one can take an 
exponentially growing sequence, e.g. pk — 2 k . This holds also in our case, in fact 
we can prove that 11.7fl is equivalent to: 




limsup > -7^ - s7V21og2 ) < +00. 

n^+oo \ - 2 

A proof of this claim can be found in |Gal| . 

When n = 2, under the above condition (non-resonance and Siegel domain), 
rescaling time by —UJ2 (assuming 0J2 ^ 0), the ODE associated to the vector field 



4 See also ICMI where a similar problem for germs of diffeomorphisms of (C, 0) has been studied. 
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can be rewritten as: 

I Z\ = LUZl + h.o.t. 



(1.8) 



i 2 = -z-2 + h.o.t. 



where lu = —loi/lu2 G (M\Q) + and high order terms means C(|z|' Q ') with |a| > 2, 
namely only the ratio of the eigenvalues enters. Then the Bruno-s condition can 
be slightly weakened (see |CMp : 




(1.9) limsup | > j - slogn | < +oo , 

n — >+oo 

where k(n) is defined by < n < Qk(n)+i an d (q n )n are the denominators of the 
convergents |HW| to u>. 

We remark that in both cases the new conditions are weaker than Bruno's con- 
dition, which is recovered when s = 0. When n = 2 we prove that the set (J B s is 
PSL (2, Z)-invariant (see remark m|l . 

The main result of |Ca2| in the case of Gevrey-s classes reads: 



Theorem 1.3 (Gevrey-s linearization). Let u>i, ...,ui n be real numbers and A = 
diagiiuJi , . . . , ico„); let D\ — {z G C™ : \z{\ < 1 , 1 < i < n} be the isotropic polydisk 
of radius 1 and let F : D\ — > C™ be an analytic function, such that F(z) = Az+f(z), 
with /(0) = Df(0) = 0. If A is non-resonant and verifies a Bruno-s, s > 0, 
condition (|1.7(l ( or condition (|1.9|l if n — 2), then there exists a formal Gevrey-s 
linearization H . 

The aim of this paper is to show that the Gevrey character of the formal lin- 
earization can give information concerning the dynamics of the analytic vector field. 
Let F(z) = Az + f(z) verify hypotheses of Theorem 1 1.81 assume moreover Xp not 
to be analytically linearizable. We will show that even if there is not a Stable do- 
main, where the dynamics of Xp is conjugate to the dynamics of its linear part, 
we have an open neighborhood of the origin which "behaves as a Stable domain" 
for the flow of Xp for finite but long time, which results exponentially long: the 
effective stability GFGHI IGPj of the equilibrium solution. 

In the case of analytic linearization, |i?j(z)|, j = l,...,n, is constant along the 
orbits, namely it is a first integral and the flow of Ijl.lfl is bounded for all t and 
sufficiently small \z \. 

We will prove that any non-zero zq belonging to a polydisk of sufficiently small 
radius r > 0, can be followed up to a time T — 0(exp{r^ 1 ^ s }), being s > the 
Gevrey exponent of the formal linearization, and we can find an almost first integral: 
a function which varies by a quantity of order r during this interval of time. More 
precisely we prove the following 

Theorem 1.4. Let n G n > 2. Given realoji, . . . , uj n consider A — diagiiui , . . . , ioj r . 
let F : D\ — > C" be an analytic function, such that F(z) — Az + f(z), with 
/(0) = Df(0) — 0. If A is non-resonant and verifies a Bruno-s, s > 0, condi- 
tion (|1.7fl (or (|1.9|) if n = 2), then for all sufficiently small < r** < 1, there 
exist positive constants A**, C„* such that for all < \zq\ < r»„/2, the so- 
lutions z{t; zq) are well defined and verify \z(t; Zq)\ < C»»r»», for all \t\ < T, = 

A~} exp^B^ (r«/|zo|) 1,/s j. 
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We want to stress here that when s — > the stability time goes to infinity, 
because the critical exponent of stability time is 1/s: namely one get stability. At 
the same time the Bruno-s condition "tends" to the classical Bruno condition, 
which is a sufficient condition to ensure analytic linearizability and hence stability 
under our assumptions. So we recover the classical stability result as limit of longer 
and longer effective stability times. 

Results similar to Theorem 11.41 have been obtained in |GFGS| for hamiltonian 
vector fields and in |GP| for real reversible systems of coupled harmonic oscillators. 
In both papers effective stability is proved by assuming the linear part of the vector 
field to verify some Diophantine condition 5 CD (7, r), for some 7 > and r > n— 1, 
and the critical exponent of stability time is 1/r. In our result, too, the critical 
exponent of stability time depends on some arithmetical property of the linear part 
of the vector field but in a more general way in fact we assume A 6 B s D CD (7, r), 
for all 7 > and r > n — 1. 

The second remark is that in GF(?51 IGP| effective stability is obtained using 
some partial normal form, then working on it and using the Poincare summation 
at the smallest term (see Lemma T2.2JI . the proof is done. Here the method used is 
completely different: we first linearize formally the system and then using properties 
of the formal linearization we conclude still using the Poincare summation at the 
smallest term. This method introduce also our main drawback: we must assume A 
to be non-resonant (to linearize) and this prevents us from considering real vector 
fields and hamiltonian ones, where an "intrinsic" resonance is present. 

In section we discuss the relation between the Bruno-s condition and other 
arithmetical conditions. 



In this part we will prove our main result, Theorem ll.4l The proof will be divided 
into three steps: first we use the Gevrey-s character of the formal linearization H, 
given by Theorem 11.31 to find an approximate solution of the conjugacy equa- 
tion (|1.3|) up to a (exponentially) small correction (paragraph 12. 1|) : then we prove 
a Lemma allowing us to control how the small error introduced in the solution 
propagates (paragraph s. 2J1 . Finally we collect all the informations to conclude the 
proof fparagraph !2.3|l . 

2.1. Determination of an approximate solution. Let F verifies hypotheses of 
Theorem II .41 and let us consider the first order differential equation in C™, n > 2: 



By Theorem II .31 this system can be put in linear form by a formal power series H 
which belongs to C s and it solves (formally): 



we observe that one can choose H{z) = z + 0(\z\ 2 ). 



A vector lj £ K n belongs to CD(-y,r) if there exists 7 > and r > n — 1 such that for all 
a £ N™ and all j £ {1, . . . , n} one has: \a ■ ui — uij\ > 7|a| — T . Let A = diag{u)i , . . . , ui n ), then A 
verifies a Diophantine condition if lu does. 



2. Proof of the main Theorem 



(2.1) 




(2.2) 



±H(z)=AH(z), 
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Since H = h a z a G C s , there exist positive constants A\ and B\ such that 

(2.3) \h a \ <^sr |Q| (l«|!) s V|«|>1. 

For any positive integer iV we consider the vectorial polynomial, sum of homoge- 
neous vector monomials of degree 1 < I < N, defined by: Hn(z) — Yli=x Yl\a\=l h a z a 
and the Remainder Function: 

(2.4) K N (z) = DH N {z) ■ F(z) - AH N (z) . 
Clearly TLm{z) doesn't solve the linearization problem, but: 

(2.5) j t H N (z) = AH N (z)+K N (z), 

hence the remainder function gives the difference from the true solution and the 
approximate one. 

The following Proposition collects some properties of the remainder function. 

Proposition 2.1. Let TZn(z) be the remainder function defined in (|2.4|l and let 

a E N n , then: 

1) d?K N (0) = if\a\ < N. 

2) For all < r < 1 there exist positive constants A2 and B2 such that if 
\a\>N+ 1, then: 

^rdfn N (0) <A 2 r-^B 2 - sN (N\) s . 

3) For all < r < 1 and \z\ < r/A there exist positive constants A^,B^ such 
that: 

/I I x N+1 

(2.6) \K N (z)\ < A 3 Bz sN (N\y ' 



Where we used the compact notation jjd° 



r 
1 



ai!...Q„! a" 1 . ..a?- 



Proof. Statement 1) is an immediate consequence of the definition of 1Z n- 

To prove 2) we observe that TZn(z) is an analytic function on D\, being obtained 
with product of analytic functions, then one gets by Cauchy's estimates for all 
< r < 1 and for all \a\ > N + 1: 



f 



1 1 



(2-7) -fits® < (27r) „ rM+1 mg\DHs-F{z)\, 

because d a HN = for \a\ > N + 1. Recalling the Gevrey estimate (|2.3|l for Hn 
and the analyticity of F we obtain: 



(2.8) ^-d^TZ N {0) <A 2 B2 sN (N\) s r-W 



for some positive constants A2 and B2 depending on the previous constants, on the 
dimension n and on F. 

To prove 3) let us write the Taylor series 1Zn{z) — Y^\ a \>N+i 't\^z'^N{^)z a : the 
bound on derivatives (|2.8() implies the estimate (|2.6|l for all \z\ < r/4 and for some 
positive constants A3 and B3. □ 

The bound 1)2. 6JI on TZn(z) depends on the positive integer N, so we can deter- 
mine the value of N for which the right hand side of 1)2. 6[) attains its minimum, 
that's Poincare's idea of summation at the smallest term. 
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Lemma 2.2 (Summation at the smallest term). Let TZn(z) defined as before and let 
< r* < 1/4 then there exist positive constants A&, B4 such that for all < \z\ < r* 
we have: 

(2.9) \nx(z)\<A A exp{-Bj^-\ }, 

where N — [B^ (f**/|z|) J and \x\ denotes the integer part of x G R. 

Proof. Let us fix < r* < 1/4, then for < \z\ < r* by Stirling formula we obtain: 

(2.10) 1^(^)1 < A 4 (A^i?^ 1 (k|/r,) 1/s ) e~ sN , 

for some positive constant A^. The right hand side of (|2.10l) attains its minimum 
at N — B% [r / \z\) x ^ s 1 evaluating the value of this minimum we get (|2.9() with 
B± = B 3 . □ 

2.2. Control of the "errors". Let us define H(z) = H^(z) and IZ(z) — TZ^(z), 
being N the "optimal value" obtained in Lemma 1^21 We remark that H(z) doesn't 
solve 112.211 but the "error", TZ(z), can be made very small: exponentially small. We 
will prove that for initial conditions in a sufficiently small disk, one can follow the 
flows for an exponential long time without leaving a disk comparable size. 

Lemma 2.3 (Control the flow). Let a,b,a and R be positive real numbers. Let 
T = Ra~ 1 e b ^ 2 ^ and let us consider the Cauchy problem: 

\x(0)=R . 
Then < x(t) < 2R for all \t\ < T. 

Proof. Let us write the Cauchy problem in integral form: 



(t) = R+ f ae- b/ ( x{s))a ds, 
Jo 



x(t) is trivially monotonically increasing, hence the same holds for the function 
t \— > ae~ b /( x ( t ^ a . Let us suppose that there exists < to < T, for which x(to) = 2i?; 
then 

2R = x(t Q ) =R + [° ae- b/(x(s)]a ds<R + t ae- b/{x{ta))a , 



namely t > Ra 1 e fc /( 2i? '" = T, which gives a contradiction. Hence either x(to) > 
2R for all < t < T or x(t ) < 2R, but the first case have to be excluded because 
x(0) = R < 2R. 

The case t < can be handle in a similar way by showing that 1 1— > x(t) doesn't 
decrease too much. □ 

Let r* as in Lemma l2~2l define p(z) = \H.(z)\ for all < \z\ < r*, then Lemma l2~2l 
admits the following Corollary, which allows us to control the evolution of p(z). 

Corollary 2.4. Let < r* < 1/4, then there exists < r*„ < r» and positive 
constants A^^B* such that for all < \z\ < r*„ we have: 

(2.11) l^)^ A M- B *{j(^j ,s }- 
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Proof. Let j £ {1, . . . ,n} and let us consider the time evolution of \Hj{z(t))\. If 
7i was a solution this would be a constant of motion, this is not the case but its 
evolution is nevertheless very slow. In fact thanks to (|2.5|l we get: 

^■(z)|<|7^(z)|, 

hence a similar statement holds for p{z) — \TL{z)\ — maxKjXn \TCj(z)\. 

We want now to express the exponential smallncss of |7?.(z)| in terms of p(z) in- 
stead of \z\. H(z) is tangent to the identity close to zero and then locally invertiblc. 
The inverse is still tangent to the identity, vanishing at zero and analytic in a neigh- 
borhood of the origin, then sufficiently close to the origin we have \TC^ 1 (w)\ < C\uu\, 
for some C > 0. Finally we can take \z\ sufficiently small, say \z\ < r** for some 
< r** < r*, s.t. \z\ < C\H(z)\ < r*, hence: 

\K{z)\ < A 4 exp{-B 4 (j^-j } < A 4 exp{ -B A C~ xl 

and the claim follows with A* = A 4 and £?» = B 4 C~ X I S . □ 

2.3. End of the proof. We are now able to conclude the proof of the main The- 
orem II .41 Take any < \z \ < and let po = |W(z(0; z o))\, then there exists a 
positive constant C\ s.t. po < Ci\zq\. By Corollarv l2.4l we have 

(2.12) J t P^ ^ A * cx p{ - B * ( r */p( z & z o))f S } ■ 

Let us call R = po, a = A*, b = B*rl^ s and a = l/s then we can apply Lemma l2.3l 
to conclude: 

(2.13) p(z(t; z )) < 2p < dr,* V|t| < T» = poA' 1 exp^B* 

This implies that TL(z{t] zq)) is well defined in this interval of time, it is not constant 
and \Tt(z(t; z ))\ << 3Cir M . Recalling that H(z) is tangent to the identity close 
to zero, we have \z\ < C 3 \H(z)\ for some positive C 3 . 

Then setting A** = 2A*r~J-, B** = B* (r„/(2r*»)) 1 ^ and C** = C1C3, we get: 

\ z (t'i z o)\ — C**f** , 
for all \t\ < A-^exp[B„ (j^)'' }■ 

3. Arithmetical conditions 

In this paper we proved that any analytic germs of vector field of (C",0) with 
diagonal, non-resonant linear part has an effective stability domain, i.e. stable up 
to finite but "long times", close to the stationary point, provided the linear part 
verifies an arithmetical Bruno-like condition depending on a parameter s > 0, 
which in the case of 2 dimensional vector fields can be put in the form: 

/fc(n) \ 
t / log <? 7 + l . \ 

limsup y, slog n < +00 . 
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Remark 3.1 (Invariance of |J S >0 B s , n = 1 under the action of PSL(2, Z))- The 
continued fraction development jHWI IMMY) of an irrational number uo gives us the 
sequences: (a k ) k >o and (uj k ) k >o- Then we introduce (0k)k>-i defined by f3-\ = 1 
and for all integer k > 0: f3 k — Wj—Q^ki which verifies : 1/2 < PkQk+i < 1 and 
QnPn-i + Qn-iPn = 1, where q k 's are the denominators of the continued fraction 
development of to. We claim that condition Bruno-s (|1.9f) is equivalent to the 
following one: 

k 



(3.1) limsup /3j-i \ogLUj 1 + s\ogf3 k -i < 



k— »-j-oo 



This can be proved by using the relations between Pi and qi, to obtain the bound, 
for all integer k > 0: 



£ (fh-i logo* + I < £ |A-i log Aa+x 



a? 

;=o v ^ ' 2=0 



A— i logA- 



Qi — 1 o , 



< 18. 



where we used the convergence of series ^2q t 1 and ^2q t 1 log a; f see MMY page 

To prove the invariance of \J S B s under the action of PSL(2, Z), is enough to 
consider its generators: Tlu = u) + 1 and Suj = 1/lu. For any irrational u), T acts 
trivially being f3 k (Tu>) = /3k(u>) for all k, whereas for S we have = L>f3k-i(Su>) 

for all k > 1. Let uj be an irrational and let u>' = oj^ 1 , let us also denote with a I 
quantities given by the continued fraction algorithm applied to uj' , then using (|3.1|l 
one can prove: 

(k \ fc+l 

X]$-i lc, g<4 1 + SUJ o ll °sP'k-i =C{uj,s)+^2/3 :j - 1 \oguj j - 1 + s\ogl3 k , 
3=0 J j=0 

where C(u>, s) = ujq (loguj^ 1 — s logt^o) + Si=o A-i ^oguf 1 , from which the claim 
follows. 

Let us consider a slightly stronger version of the Bruno-s condition: uj G (0, 1)\Q 
belongs to B s if: 

(3.2) lim fv^±l- s logJ <+oo, 

where (q n ) n are the convergents to u>, and let us introduce a second arithmetical 
condition denoted by B' s to be the set of irrational numbers whose convergents 
verify: 

(3.3) lim ^ML = 8 . 

k^+co q k \ogq k 

We state without proof the following proposition, and we refer to |Calj , to all 
details: 

Proposition 3.2. Let s > and let ui G (0,1) n B s . Then if uj is not a Bruno 
number then u> G B' s , otherwise u) G B' . 

Therefore if u) G B s \ B then the denominators of the convergent to u) can grow 
like a factorial, more precisely, qk+i = Ol (q k *) s ) , is allowed. 
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